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We investigate the three-dimensional, time-reversal invariant topological superconductors with 
generic interaction by their response to external fields. The first description is a gravitational 
topological field theory, which gives a Z2 classification of topological superconductors, and predicts 
a half-quantized thermal Hall effect on the surface. The second description introduces an s-wave 
proximity pairing field on the surface, and the associated topological defects give an integer Z 
classification of the topological superconductors. 
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INTRODUCTION 



Recently, topological insulators (TI) and topological 
superconductors(TSC) have attracted great attention in 
condensed matter physics 0-0]. Topological insulators 
are states that are incompressible in the bulk but have 
topological surface states protected by the bulk topologi- 
cal properties which remain robustly metallic upon arbi- 
trary perturbation on the surface, within the given sym- 
metry class. Time-reversal invariant topological insula- 
tors in two and three dimensions (2D and 3D) have been 
theoretically predicted and experimentally realized in re- 
cent vears(5Hll|. Topological superconductors are ana- 



log of topological insulators in superconductors, which 
have a full superconducting gap, and gapless Majorana 
edge states propagating on the boundary. Time-reversal 
breaking topological superconductors in 2D 12j and time- 
reversal invariant topological superconductors in 2D and 
3D (l3rfl5| have been theoretically predicted but have not 
been verified experimentally. The most promising candi- 
date for 3D topological superconductor (more precisely, 
superfluid) is the 3 He-B phase. 

There are two approaches for the classification and 
description of TIs and TSCs: topological field the- 
ory (TFT) and topological band theory (TBT). In the 
TFT approach, topological insulators are defined by the 
quantized responses derived from low-energy effective 
actions 16]. TFT also reveals the relations among TIs in 
various dimensions. For instance, the 4D Chern-Simons 
insulator [17} is the root state of 3D and 2D time re- 
versal invariant (TRI) TIs. The 3D topological insula- 
tor is characterized by the topological effective theory 
S = ^- F ^ F or^ vaT with = 7T mod 2tt, which de- 
scribes the electromagnetic response of the topological 
insulator when time-reversal symmetry is broken on the 
surface. [H El The TBT was first developed in the work 
of Kane and Mele 3, 11 on the Z2 topological classifi- 
cation of quantum spin Hall insulator, which was later 
generalized to 3D TRI TIs 21-21 an d more recently to 



all ten symmetry classes [1J, [24[. It was also shown that 
TFT reduces to TBT in the non-interacting limit for 
time-reversal invariant TIsfioT]. 

Interestingly, in ID it has been shown that the 
non- interac ting classifications are sometimes unstable to 
interaction. 26]. Therefore, in order to define TSCs in 
dimensions higher than 1 , it is helpful to apply the topo- 
logical field theory (TFT) approach of TIs[16l Il7j . since 
the TFT describes physical response properties of the 
system, which is applicable to interacting systems. All 
topological insulators can be described in this way by 
TFT[ll, including the 2D and 3D TIs. An explicit ex- 
pression of the topological invariant of these classes in 
interacting systems are obtained in the form of Green's 
function 27]. However, there are some symmetry classes 
in which this approach does not apply. For example, the 
3D time-reversal invariant TSC is one of them. Due to 
the absence of charge conservation symmetry, no TFT 
can be obtained by the coupling to the electromagnetic 
field. 

In this paper, we introduce two new approaches to 
the classification of 3D interacting time-reversal invari- 
ant TSCs. The general idea behind that can generate a 
physical observable topological response of the system. 
Such topological response can then be used to classify 
the topological states. TIs and TSCs in different topo- 
logical classes and different dimensions require different 
external fields to probe their topological response proper- 
ties. We present two different approaches with different 
probe fields used, which applly to different physical set- 
tings. The first approach is the gravitational TFT, which 
describes the topological response of the TSC to a grav- 
itational field. Since the coupling to gravitational field 
does not require charge conservation, the gravitational 
TFT can be applied to TSC similar to the electromag- 
netic TFT approach of 3D Tl[l6j]. The topological term 
is the gravitational analogy of the #-term for the gauge 
field, which is the Pontryagin invariant of the spacetime 
manifold [2I. However, similar to the electromagnetic 9 
term for 3D TI, the gravitational TFT only provides a 
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Zi classification of the 3D TSC. In the non-interacting 
limit, the nontrivial class with 9 — -k corresponds to the 
TSC states with odd topological quantum number, while 
the trivial class with 9 = corresponds to those with 
even topological quantum number. 

A natural question is whether this indicates that with 
interaction only a Z^ classification is stable, or there are 
other ways we can define the Z classification in interact- 
ing TSC. To answer this question we propose the second 
approach. We consider the proximity effect of the surface 
of 3D TSC with an s-wave superconductor. The proxim- 
ity effect introduce an s-wave pairing field to the surface 
states. When the pairing field has a 7r phase domain wall 
on the surface, there are N number of chiral Majorana 
fermions propagating on the domain wall. Since the chi- 
ral Majorana fermions are stable without any symmetry 
requirement, and can be characterized by the thermal 
current it carries, one can use the number of Majorana 
zero mode as a generic definition of TSC in 3D. From this 
approach we conclude that the integer classification of 3D 
TSC remains robust when interaction is considered. This 
approach also directly leads to an experimental proposal 
for the measurement of the topological quantum number. 
Moreover, the idea of defining interacting TIs and TSCs 
by studying the defects obtained from some symmetry- 
breaking field configuration can be generalized to generic 
symmetry classes. 



GRAVITATIONAL TOPOLOGICAL FIELD 
THEORY OF 3D TSC 

First we consider gravitational TFT of time-reversal 
invariant (TRJ) TSC in 3D. For simplicity, we first inves- 
tigate continuous models using the example of the 3 He-B 
phase as a 3D TRI TSC. 

In the fl at sp ace-time, the Hamiltonian of 3 He-B phase 
is given by[13j, [29j 



H(p) = m(p)r 3 + Act • pr 1 



(1) 



in the Nambu space ip = (cf, c±, ct , — cl), where f, 1 de- 
note spin orientation, a 1 ' 2,3 and r 1,2,3 labels the Pauli 
matrices in spin and particle-hole spaces, respectively. 
m(p) = p 2 /2M — /i is the kinetic energy term. The weak 
pairing phase with M > 0, fJ, > is topologically non- 
trivial and corresponds to the physical 3 He-B phase [13? 

while the strong pairing phase with fi < is a trivial 
superfluid phase. After rescaling the momenta Api — > pi, 
the corresponding Lagrangian can be written as 



ip[p^ + m{pi)]ip 



(2) 



where \i = 0, l,2,3(We use (0,1,2,3) interchangeably 
with (t, x, y, z)), r° = r 3 and T l = ia l r 2 ;i = 1, 2, 3. This 
is exactly the Majorana-Dirac Lagrangian if m(pi) = m 
is independent of pi . 



Now we turn to Lagrangian of 3 He-B in curved space- 
time. It is convenient to take the space-time as a closed 
4-manifolds M. To consistently define the Lagrangian 
globally, M is required to be a spin manifold [13] . Locally 
the low energy action is 



Sib = 



M 



c^fp[zT a e^ + + n# (3) 



where S Q( g = |«[F Q , is the generator of Lorentz trans- 
formation in the spinor representation, and is the 
vielbeins[28j]. Physically we can interpret the gravita- 
tional field as an external source coupling to the energy- 
momentum tensor of the fluid. In a deeply insightful pa- 
per, Volovik points out that the order parameter of the 
3 He-B phase also couples to fermions like a gravitational 
vielbein[3l|. Therefore, we can also interpret as the 
internal order parameter of the 3 He-B phase. 

Because the Majorana fermions are gapped, we can 
integrate out them to obtain a gravitational effective ac- 
tion. The term we are interested in is the gravitational 
theta term, which can be easily obtained by calculat- 
ing chiral anomaly due to gravity. Due to the non- 
invariance of fermion Jacobian, a chiral transformation 
ip — > exp(ir 5 #/2)-!/>, which inverts m to —m, gener- 
ates the topological term in the effective action 32, |33| 



as So = 7}f M A(M), where A(M) is the Dirac genus 



of M[28|. Since the Hamiltonian with positive mass 
m = — fi > is adiabatically connected to a trivial su- 
perconductor by taking fi — > — oo, by a chiral rotation 
of 9 = 7r we can obtain the action of the 3 He-B phase 
with fi > as Sg = § f M A(M). It is worth noting 
that Majorana fermion has a 1/2 coefficient compared to 
Dirac fermion. More explicitly, we have the Pontryagin 
invariant 
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(4) 



with 9 = it. For a generic 4- manifold, Sg = n9/A8 with 
integer n. However, for spin manifolds the value of Dirac 
genus is restricted to J M A(M) = 2n with integer n[34j|. 
so that S$ = n9, as is expected from time reversal invari- 
ance. On a domain wall between 9 = and 9 = tt, we 
get the Chern-Simons term from the theta term as 

If 2 

S 2 d = --7- / d 3 xe^ p Tr(uj^d v uj p + -ui^,ui u cj p ) (5) 

which describes the surface of 3 He-B phase in the pres- 
ence of surface magnetization. 

For a generic TSC with a non-relativistic Hamilto- 
nian, the direct coupling to gravity may be complicated. 
However, physically one can always consider the non- 
rclativistic Hamiltonian as a low energy limit of an under- 
lying relativistic theory, so that the coupling to gravity is 
in principle always well-defined, and the topological term 
has a coefficient of 9 = or it mod 2ir as is required by 
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TR symmetry. The topological term ((4]) is insensitive to 
an arbitrary deformation of the fermion Hamiltonian as 
long as the fermions remain gapped and TRI is preserved. 
For a non-interacting TSC with topological number N, 
its Hamiltonian can always be deformed to N copies of 
the 3 He-B Hamiltonian 2J| , so that the topological term 
has = Nir. However, for even integer N, exp(iSe) = 1 
and Sg has no physical effect. We conclude from this fact 
that the gravitational TFT gives a Zi classification of 3D 
TRI TSC, which is weaker than the integer classification 
of the non-interacting system. Because gravitational re- 
sponses are related to thermal responses [35], the gravi- 
tational TFT predicts topological thermal responses, as 
we shall discuss in the next section. 



SURFACE STATE DESCRIPTION 

To understand the difference between Zi and Z clas- 
sification, it is helpful to study the topological surface 
states. For the 3 He-B phase, the surface state is a single 
flavor of 2D massless Majorana fermion 13, [iij]. For our 
discussion here, it is most convenient to use the Majorana 
basis in which the surface state effective Hamiltonian is 
given by 



(6) 



with 77k = ?}_k the two-component Majorana fermion 
field. The two components carry opposite spin, and the 
time-revesal operation is defined T _1 ?7kT = io~ y 7]-k. It 
should be noticed that only Pauli matrices a z and o~ x 
are allowed in this Hamiltonian since the single particle 
Hamiltonian /i(k) = v (a z k x + a x k y ) satisfies the anti- 
symmetry condition h{— k) = — /i T (k). No mass term is 
allowed by time-reversal symmetry since the only possi- 
ble mass term io~ y breaks time-reversal symmetry. More 
generically for a topological superconductor with topo- 
logical quantum number TV, the surface state consists of 
multiple copies of the Majorana fermions with different 
chirality: 



H = ^ r]J_ k v (a 7 k x + a x k y ) ri s ^ 

s=l k 

N- 

+ X!X!^-k t; ( a z k * ~ °xky) Ips.k (7) 
s=l k 

where r] s ^ and i/j s ,k are the "left-handed" and "right- 
handed" Majorana fermions, and the integer N + — 7V_ = 
N is determined by the bulk topological invariant. Such a 
definition of "chirality" for Majorana fermion is possible 
because a y is not allowed to appear in the linear- k terms, 
therefore the two kinds of Majorana fermions 77k and "0k 
cannot be deformed to each other. 



More rigorously, the topological invariant N = 
N- can be defined by the following formula: 



iV=-Wlnd (r 



dh(k) dh(k) 

dk x 



dk v 



(8) 



k=K ; , 



with T = idyl the time-reversal transformation matrix 
for all Majorana fermions, and K; runs over all time- 
reversal invariant momenta (0, 0), (0, 7r), (-tt, 0), (-7T, n) in 
the surface Brillouin zone. Ind(Af) for a Hermitian ma- 
trix M is defined as the number of positive eigenvalues 
minus the number of negative eigenvalues. 

From the discussion above we have seen from the sur- 
face state point of view why the topological classifica- 
tion is Z for non-interacting theory. The relation of 
this surface state picture to the gravitational TFT (|4| 
is similar to that between the Dirac surface state and the 
topological electromagnetic action of 3D Tl[l6j]. Con- 
sider a generic T-breaking mass term of the surface states 
H m = J2 S m sVj^yVsk + J2 S m's^J-kVy^sU- With the 
mass terms the surface states are gapped and we obtain 
a quantized thermal Hall coefficient 



Trk 2 B T 
24h 



'sgn(TO s ) 



E 

s=l 



sgn(m^ 



(9) 



Each Majorana cone contributes half the thermal Hall 
coefficient of that of a 2D (p + ip) superconductor [1 2l 136 - 
38] . Such a thermal Hall effect on the surface is consis- 
tent with the description of the bulk TFT g]) since the 
Pontryagin invariant can be reduced to a gravitational 
Chcrn-Simons term on the boundary as eq.([5]). However, 
due to the dependence on the sign of mass terms, the 
thermal Hall coefficient is only determined by the sur- 
face state theory ((7]) mod 2 x 1 This is consistent 
with the observation that the coefficient of the bulk TFT 
(U]) takes Z2 value or tt mod 2ir. In other words, only 
the topological superconductors with odd N necessar- 
ily have thermal Hall effect on the surface, while those 
with even N can have zero thermal Hall coefficient for 
some choice of mass terms. It's interesting to note that 
the same gravitational TFT approach can be applied to 
time-reversal invariant TI, for which the surface thermal 
Hall conductivity will be doubled. Thus we see that the 
gravitational TFT description is trivial for 3D TI. 



INTEGER Z CLASSIFICATION OF 3D TRI 
SUPERCONDUCTOR BY TOPOLOGICAL 
DEFECTS 

The analysis above seems to indicate that only the 
topological superconductors with odd topological quan- 
tum number are stable. However, in the following we 
will show that the Z classification is actually robust 
with interaction, although it cannot be captured by the 
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gravitational TFT. As has been discussed in the in- 
troduction part, a general philosophy of classifying in- 
teracting topological phases is by defining a physical 
topological response of the system to a proper exter- 
nal field. The gravitational TFT can only probe the 
Z2 part of the topological invariant, but the full Z can 
be probed by another physical field — an s-wave pair- 
ing field. We consider the generic 3D time-reversal in- 
variant Majorana fermion system with a Hamiltonian 
H = J d 3 rd 3 r'ri T (r)h(r, r')ry(r') + H int . Here h is an 
anti-symmetric Hermitian single-particle operator, and 
Hi nt describes non-quadratic terms in Majorana fermion 
r\. The Hamiltonian is time-reversal invariant satisfying 
T~ 1 HT = H, with the anti-unitary time-reversal opera- 
tor defined by T~ 1 ry(r)T = Tr](r) and T 2 = — 1. Now we 
introduce the following perturbation term to the Hamil- 
tonian 

H A = -if d 3 rA(r)?j r (r) [T _1 7y(r)T] 

= J d 3 riA(r))j T (r)rr/(r) (10) 

with A(r) being a real field. To see the physical mean- 
ing of this term, we can take the 3 He-B Hamiltonian (JTJ 
as an example. In the complex fermion basis = 
— J d 3 riA(r) (cj k (r)cj k (r) — h.c.^j which is an s-wave 
pairing term with imaginary pairing order parameter. It 
is important to notice that such a T-breaking pairing 
term pop is completely determined by the representa- 
tion of time-reversal transformation and is independent 
of the basis choice and any other detail of the system. 

Now we consider the time-reversal invariant TSC with 
the perturbation (fTU|) only near the surface. Physically, 
such a field can be induced by the proximity effect be- 
tween the TSC and an s-wave superconductor on top of it. 
Since this term breaks time-reversal symmetry, the sur- 
face states of the TSC can become gapped. The effect of 
this s-wave pairing field to the surface state can be seen 
straightforwardly from the surface effective theory ([7]). 
Since T~ l, q s T = iu y r] Sl T~ x i\)tT = ia y ipt for the surface 
Majorana fermions, the s-wave pairing term (|10j) is re- 
duced to HAsurf = J cPrA(r) (J2 S vl^yVs + J2t ^t^y^t) 
on the surface. In other words, an uniform pairing field 
A(r) = A > corresponds to a uniform mass term 
m s = m' t = A for all the surface Majorana fermions. 
Consequently the thermal Hall conductivity given by Eq. 
([9]) is K xy — N w ith N the bulk topological invari- 
ant, so that by introducing the pairing field, the topolog- 
ical quantum number iV can be obtained without ambi- 
guity. 

The statement above can be made more explicit by 
considering a domain wall on the surface where A(r) 
changes sign, as shown in Fig. [TJ Physically, such 
a domain wall is a 7r-Josephson junction of the s-wave 
superconductor on top of the time-reversal invariant 




FIG. 1: Illustration of the experimental setting for the mea- 
surement of the integer-valued topological quantum number 
of 3D TSC. s-wave pairing is introduced on the surface by 
proximity effect, with a 7r-Josephson junction. The topolog- 
ical quantum number N of the bulk is characterized by the 
number of chiral Majorana fermions propagating along the 
junction illustrated by the red arrows, which carry a quantized 

thermal current I E — N ^ 2^ f° r l° w temperature T — > 0. 

TSC. For the surface Majorana fermions, the mass term 
changes sign across the domain wall, therefore there are 
chiral Majorana fermion propagating along the domain 
wall [l2j . On comparison, similar 7r-Josephson junction 
on the surface of topological insulator leads to non-chiral 
Majorana fermions [39]. Due to the uniform mass for 
all Majorana fermions, the number of chiral Majorana 
fermion is given by iV+ while the number of anti-chiral 
Majorana fermion is given by iV_. Consequently there 
are N = N + — N- copies of chiral Majorana fermions 
which are robust gapless states on the domain wall. It is 
important to notice that a chiral Majorana fermion in ID 
is well-defined and robust even in an interacting system, 
since it has a chiral central charge c — 1/2 and thus car- 
ries a chiral thermal current I E = Thus one can 
use the chiral thermal current as a physical observable 
way to define chiral Majorana fermions, and the chiral 
Majorana fermions on the domain wall can be used as a 
generic definition of the bulk 3D TSC as follows: 

• For a generic 3D time-reversal invariant supercon- 
ductor, apply a small pairing field in the form of 
Eq. (|10[) with a domain wall of A(r), as shown in 
Fig. [TJ The superconductor is a TSC with topolog- 
ical quantum number N € Z if the thermal current 
propagating along the domain wall at temperature 
T is I E = N^^- in the T -> limit. 
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